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ABSTRACT 

We study the effects of channel coupling in the excitation dynamics of 
giant resonances in relativistic heavy ions collisions. For this purpose, we 
use a semiclassical approximation to the Coupled-Channels problem and 
separate the Coulomb and the nuclear parts of the coupling into their main 
multipole components. In order to assess the importance of multi-step 
processes, we neglect the resonance widths and solve the set of coupled 
equations exactly. Finite widths are then considered. In this case, we 
handle the coupling of the ground state with the dominant Giant Dipole 
Resonance exactly and study the excitation of the remaining resonances 
within the Coupled-Channels Born Approximation. A comparison with 
recent experimental data is made. 


1 Introduction 

Relativistic Coulomb Excitation (RCE) is a well established tool to unravel interesting 
aspects of nuclear structure [Q. Examples are the studies of multiphonon resonances 
in the SIS accelerator at the GSI facility, in Darmstadt, Germany p|, Important 
properties of nuclei far from stability have also been studied with this method. 
The RCE induced by large-Z projectiles and/or targets, often yields very large cross 
sections, which depend on the nuclear response to the acting electromagnetic fields. 
In such cases, one expects a strong coupling between the excited states. This coupling 
might be responsible for the large discrepancies between experimental data of RCE 
and the calculations based on first-order perturbation theory BBa. or the harmonic 
oscillator model. 

In the present paper, we apply a semiclassical method ^ to the coupled-channels 
(CC) problem and study RCE in several collisions between heavy ions. In this method, 
the projectile-target relative motion is approximated by a classical trajectory and the 
excitation of the Giant Resonances is treated quantum mechanically |^, |^. The use 
of this method is justified due to the small wavelenghts associated with the relative 
motion. In section 2, we neglect the resonance widths and introduce the semiclas¬ 
sical CC-equations for relativistic Coulomb excitation. The time-dependent matrix- 
elements of the main multipole components of the Coulomb (section 2.1) and nuclear 
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(section 2.2) parts of the coupling interaction are calculated. The CC-equations are 
then solved (section 2.3) in some limiting cases. Section 3 is devoted to the excita¬ 
tion of resonances of hnite widths. Generalizing the schematic treatment of ref. 0, 
we present an “exact” solution for the coupling between the g.s. and the domi¬ 
nant GDR. The excitation of the weaker resonances are then evaluated through the 
Goupled-Ghannels Born Approximation (GGBA), from the g.s. and GDR amplitudes. 
In section 4 we apply the results of the previous sections to specihc cases and make 
a comparison with recent experimental data. Finally, in section 5, we summarize our 
results and present the conclusions of this work. 


2 The semiclassical method for the CC-problem 

In relativistic heavy ion collisions, the wavelength associated to the projectile-target 
separation is much smaller than the characteristic lengths of system. It is, therefore, 
a reasonable approximation to treat r as a classical variable r(f), given at each instant 
by the trajectory followed by the relative motion. At high energies, is also a good 
approximation to replace this trajectory by a straight line. The intrinsic dynamics can 
then be handled as a quantum mechanics problem with a time dependent hamiltonian. 
This treatment is discussed in full details by Alder and Winther in ref. 

The intrinsic state \'ip(t) > satishes the Schrodinger equation 

[ft + r(r«))| \m >= . (1) 

Above, h is the intrinsic hamiltonian and V is the channel-coupling interaction. 

Expanding the wave function in the set {\m >; m = 0,A^} of eigenstates of h, 
where N is the number of excited states included in the Goupled-Ghannel problem, 
we obtain 

N 

IV’(^) >= ^rnit) |m > exp ( - iEmt/h^ , (2) 

m=0 

where Em is the energy of the state \m >. Taking scalar product with each of the 
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states < n|, we get the set of coupled equations 

N 

ih an{t) = ^ <n\V\m> am(t) n = 0to A^. (3) 

m=0 

It should be remarked that the amplitudes depend also on the impact parameter b 
specifying the classical trajectory followed by the system. For the sake of keeping the 
notation simple, we do not indicate this dependence explicitly. We write, therefore, 
ttnit) instead of a„(6, t). Since the interaction V vanishes as t —>• ±cx), the amplitudes 
have as initial condition a„(t —>■ —cxo) = 5(n, 0) and they tend to constant values as 
t —> oo. Therefore, the excitation probabity of an intrinsic state |n > in a collision 
with impact parameter b is given as 

Pn{b) = \an{oo)\^ . (4) 

The total cross section for excitation of the state In > can be approximated by the 
classical expression 

cTn = 27r J Pn{b) bdb . (5) 

Since we are interested in the excitation of specihc nuclear states, with good an¬ 
gular momentum and parity quantum numbers, it is appropriate to develop the time- 
dependent coupling interaction V{t) into multipoles. In ref. pj, a multipole expansion 
of the electromagnetic excitation amplitudes in relativistic heavy ion collisions was 
carried out. This work used hrst order perturbation theory and the semiclassical 
approximation. The time-dependence of the multipole interactions was not explicitly 
given. In section 2.1 we show how this time-dependence can be explicitly obtained, 
from the Taylor-series expansion of the Lienard-Wiechert potentials [Q and the con¬ 
tinuity equation for the nuclear current. 

In section 2.2 we deduce the time-dependence and the multipole decomposition of 
the nuclear interaction in relativistic nucleus-nucleus collisions. The nuclear absorp¬ 
tion at collisions below grazing impact parameter is also accounted for. 
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2.1 Coulomb excitation 


We consider a nucleus 1 which is at rest and a relativistic nucleus 2 which moves 
along the 2 ;-axis and is excited from the initial state \IiMi > to the state \IfMf > 
by the electromagnetic held of nucleus 1. The nuclear states are specihed by the 
spin quantum numbers h, If and by the corresponding magnetic quantum numbers 
Mi and Mj, respectively. We assume that the relativistic nucleus 2 moves along a 
straight-line trajectory with impact parameter b, which is therefore also the distance 
of the closest approach between the center of mass of the two nuclei at the time t = 0. 
We shall consider the situation where b is larger than the sum of the two nuclear radii, 
such that the charge distributions of the two nuclei do not strongly overlap at any 
time. The electromagnetic held of the nucleus 2 in the reference frame of nucleus 1 is 
given by the usual Lorentz transformation [Q of the scalar potential ^(r) = Zie/|r|, 
i.e., 


A(r',t) 


7 0 

V 

- 7 


b' — b, 7(2;' — vt) 
b' — b, 7 (z' — vt) 


( 6 ) 


c L 

Here b (impact parameter) and b' are the components of the radius-vectors r and r' 
transverse to v. 

The time-dependent matrix element for electromagnetic excitation is of the form 


Vjft) =< !,M, 



V 

72 



0(r', t)\IiMi > 


( 7 ) 


A Taylor-series expansion of the Lienard-Wiechert potential around r' = 0 yields 


0(r', t) = 70[r(t)] + 7V0[r(f)] ■ r' + ■ ■ ■ 


( 8 ) 


where r = (b, ynt), and the following simplifying notation is used: 


V 0 [r] 


V' 0 (r', t) 

r '=0 


( 9 ) 
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Thus, 


Vfi{t) =< IfMfl p(r')-- ■ J(r') 70(r) + 7 r' ■ V0(r) \IiMi > . (10) 


Using the continuity equation 


V ■ J = —i cn p , 


( 11 ) 


where u = {Ef — Ei)/h, and integrating by parts, 

rV' vi r 


v,M =< i/Mf 

In spherical coordinates 


J(r) 




70(r + 'jr' ■ V0(r) 


ILM > 


r' ■ V.#. = ^ a, <■' n;. 


where 


We will use the relations 


M=-l 


cx^ C/i ■ V (j) 


r V ^ V Att 

^ ^ ^ ^ VT ' 


lo) 


and 

where L = —iV x r. 
Then, one can write 


V X = i{k + l)V{r%^) 


J • (— - [70 + 7r' ■ V0 


= -7J 


^i=-l 


^W-r') - i/^ 

a 
iu! 


( 12 ) 


(13) 


(14) 


(15) 


(16) 


X { E It V'(r'n,) V'(r' 17 „)} 


.(17) 


The last term in the above equation can be rewritten as 


(J—) ( 

r'.V0) = ^ J- ' 

V c^/ \ 

2c2 L 


eo(r'-V0) + (r'-eo)V0 




eo(r'-V0)-(r'-eo)V0 
(18) 
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The first term in this equation is symmetric under parity inversion, and contributes 
to the electric quadrupole (E2) excitation amplitudes, since 


2c2 


J 


eo(r' ■ V0) + (r' ■ eo)V(/> 


2c2 


J -V' 


z'{r' ■ V0) 


(19) 


The second term in eq. (^) is antisymmetric in J and r', and leads to magnetic dipole 
(Ml) excitations. Indeed, using eqs. (|T3|-^), one finds 


2c2 


J ■ 


eo(r' ■ V0) - (r' ■ eo)V0 


2c2 


J ■ 


Ajc ^ / X. 



^l=-l 


( 20 ) 


Thus, only the first two terms on the right-hand-side of eq. (0) contribute to the 


electric dipole (El) excitations. Inserting them into eq. (|1^), we get 


Wf'b)=7 


/ 47r _ ^ ^ 

— ^ (-1)^ (3^ < IfMf\M{El, -/r)| W, > 

/2=-l 


( 21 ) 


where 


and 


M{E 1 ,-fi) = ^ j = j <iV p(r) r y„(r) , 


( 22 ) 


d(j) 

(3± = -a/, = -(V0 ■ ^m) = 6/2 ■ ^ 

Po — —tto — ■ 


(23) 


The derivatives of the potential 0 are explicitly given by 


d(j) 

dh^. 


= Vb, 


r'=0 


= —X 


Zie 


[6^ -f 72t;2i2j3/2 


_ 2 . _ -^1^ _ 

r'=0 “ ^ [52 ^2y2^2] 3/2 


Using the results above, we get for the electric dipole potential 


(24) 


I2ir 


bf‘’(*) = \/T7 '|£i(7) MdEl,-l)-Mf,(El,l) 


V2 


7 T 


(1 + ^ 2 {t) M //(El, 0) 


(25) 


7 





















where r = 7 ^/ 6 , and 


£i{r) = 


Zie 


and £ 2 {t) = 


Zier 


(26) 


(,2 [1 ^ 2 p /2 " “ J |1 + t 2 | 3/2 

are the transverse and longitudinal eletric helds generated the relativistic nucleus with 
charge Zie, respectively. From the dehnition 


J d\ J(r).L(rr„.) , 


(27) 


and eq. ( 0 ). we find 


7"'7) = i 



27 r V 

^ c 


£i(t) 


(28) 


To obtain the electric quadrupole (E2) potential we use the third term in the 
Taylor expansion of eq. (^. Using the continuity equation, a part of this term will 
contribute to E3 and M2 excitations, which we neglect. We then find that 




(E2), 


/* 


r) = 


7 3£3(t) Mf,{E2,2) + Mf,{E2,-2) 


+ 'j 6 tSsIt) -i + Si{t) Mfi{E2,-l) + Mfi{E2,l) 

7 c 


+ ^ 67 ^ (2r^ - 1 ) £: 3 (r) - i r (1 + £:i(r) Wl/j(E2,0)| (29) 

7 c J 

where Ssi^r) is the quadrupole electric field of nucleus 1 , given by 

^3(^) = ^3 [X + r2]5/2 • 

The fields Siir) peak around r = 0, and decrease fastly within an interval At ~ 1 . 
This corresponds to a collisional time At ~ b/'yv. This means that, numerically one 
needs to integrate the Coupled-Channels equations (eq. (j^)) only in a time interval 
within a range n x At around r = 0, with n equal to a small integer number. This 
will be shown latter in connection with the calculation presented in section 4. 

Using the Wigner-Eckart theorem we can write 


Mj.(EK,A = (-1)''-"' (_((^ ) < /,M,||A<(£;A)||/.M. > 


—Mf /i Mi 


(31) 

















A phase convention for the nuclear states can be found so that the reduced ma¬ 
trix elements < IfMf\\M.{E\)\\IiMi > are real numbers [Q. For the case of giant 
resonances, sum rules are very useful to guess the values of these matrix elements. It 
is usual to use the reduced transition probability 




1 


2/j -|- 1 
1 

2/j -|- 1 


^ \<IM\M{EX, fi)\IfMf> 

MiMf 

<h\\M{EX)\\If> \ 


(32) 


in terms of which the energy-weighted sum-rules yield, for the El and E2 excitations. 


= (^) 4^ 


9 r NZ 


dvr 2 mjq AEx 


(33) 


and 

B{E2- k 


(2/4+ l) 


r 15/?2 


2// -|- 1/ mjv dvrEa; 


e^x ■ 


Z"^ jA, for isoscalar excitations; 
NZ/A^ for isovector excitations; 

( 34 ) 


where iV, Z, and A are the neutron, charge, and mass number of the excited nucleus, 
respectively. In these equations it was assumed that an isolated state with energy E^ 
exhausts the sum-rule. 

The matrix elements for the transitions between multiphonon states can be deter¬ 
mined by using the Wigner-Eckart theorem and the reduced matrix elements inferred 
from sum rules, as described in sections 2.1 and 2.2. In the case of perfect phonons, 
i.e., eigenstate solutions of the harmonic oscillator, the following relation holds for 

m-. 


the reduced matrix elements for the transition 0^1 and n — 1 —>■ n 


< n — 1||I4;/Ar,i||u > p — n I < 011111 > |^ . (35) 


The factor n on the r.h.s. is the boson enhacement factor. 


2 . 1.1 Approximate solutions 

In most cases, the first-order perturbation theory is a good approximation to calculate 
the amplitudes for relativistic Coulomb excitation. It amounts to using = 6ko on 
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the right hand side of eq. (^. The time integrals can be evalnated analytically for 
the VEiif) pertnrbations, given by eqs. (pS]), (PBD, and (|2^). The resnlt is 


“1st 





Svr Zie 
3 hvb 






7 


(36) 


where Ki {K 2 ) is the modihed Bessel fnnction of hrst (second) degree, and ^ = oifo/yn. 
For the E2 and Ml mnltipolarities, we obtain respectively, 

aSf> ^ 2 ^ ^ 


30 jhvb^ ^ [MME2,2) + MME2,-2) 

+ *7 ( 2 -^) ^i(0 [Mf^iE2,-l)+MfiiE2,l) 


VgMO -M/.(E2,0)| , 


(37) 


and 


JMl) 

“1st 


Zie 

V 3 hch 






(38) 


These expressions are the same as those obtained from the formnlae dednced in 
ref. |@]. We note that the mnltipole decomposition developed by those anthors is 
accomplished by a different approach, i.e., using recurrence relations for the Gegen- 
bauer polynomials, after the integral on time is performed. Therefore, the above 
results present a good check for the time-dependence of the multipole helds deduced 
here. 

A simplihed model, often used in connection with multiphonon excitations, is the 
harmonic vibrator model. In this model, the resonance widths are neglected and the 
Coupled-Channel equations can be solved exactly, in terms of the hrst-order excitation 
amplitudes |^. The excitation amplitude of the n-th harmonic oscillator state, for 
any time t, is given by 


“il(^) = 


“isl (^) 


a/uT 


exp 


{ - |ai„(i)|V2} , 


(39) 


where aistit) is the excitation amplitude for the 0 {g.s.) —> 1 (one phonon) calculated 
with the hrst-order perturbation theory. 
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For the excitation of giant resonances, n can be identified with the state corre¬ 
sponding to a multiple n of the single giant resonance state. This procedure has been 
often used in order to calculate the cross sections for the excitation of multiphonon gi¬ 
ant resonances. Since this result is exact in the harmonic vibrator model, it accounts 
for all coupling between the states. However, this result can be applied to studies of 
giant resonance excitation only if the same class of multipole states is involved, he., if 
one considers only electric dipole excitations, and use the harmonic oscillator model, 
one can calculate the excitation probabilities, and cross sections, of the GDR, double- 
GDR, triple-GDR, etc. Eq. ( PP|) is not valid if the excitation of other multipolarities 
are involved, e.g., if the excitation of dipole states and quadrupole states are treated 


simultaneously. In ref. a hybrid harmonic oscillator model has been used. In 
this work, it is assumed that the difference between the amplitudes obtained with 
the harmonic oscillator model and with n-th order perturbation theory is due to the 
appearence of the exponential term on the r.h.s. of eq. (PU[). This exponential takes 
care of the decrease in the occupation amplitude of the ground state as a function of 
time. As argued in ref. [Q, the presence of other multipole states, e.g., of quadrupole 
states, together with dipole states, may be accounted for by adding the hrst order 
excitation amplitudes for the quadrupole states to the exponent in eq. (|3^. This 
would correct for the flux from the ground state to the quadrupole states. In other 


words, eq. (|3^ should be corrected to read 

-in 

t) = ^ —L exp{ - ^ |ai,t(7r'A', t)^/2^ 


(40) 


The harmonic oscillator model is not in complete agreement with the experimental 
hndings. The double-GDR and double-GQR states do not have exactly twice the 
energy of the respective GDR and GQR states p, P] . Apparently, the matrix elements 
for the transition from the GDR (GQR) to the double-GDR (double-GQR) state 
does not follow the boson-rule (see end of section 3). This is borne out by the 
discrepancy between the experimental cross sections for the excitation of the double- 
GDR and the double-GQR with the perturbation theory, and with the harmonic 
oscillator model 0,13] . Thus, a Goupled-Ghannels calculation is useful to determine 
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which matrix elements for the transitions among the giant resonance states reproduce 
the experimental data. 

2.2 Nuclear excitation and strong absorption 


In peripheral collisions the nuclear interaction between the ions can also induce exci¬ 
tations. This can be easily calculated in a vibrational model. The amplitude for the 
excitation of a vibrational mode by the nuclear interaction in relativistic heavy ion 
collisions can be obtained assuming that a residual interaction U between the pro¬ 
jectile and the target exists, and that it is weak. According to the Bohr-Mottelson 
particle-vibrator coupling model, the matrix element for the transition i —/ is 
given by 


yf ""’(r) =< IsMj\U\hMi >= < IsMs\Y,,\hM. > VV(r) U,{r) (41) 

where 5\ = (5\R is the vibrational amplitude, or deformation length, R is the nuclear 
radius, and U\{r) is the transition potential. 

The deformation length 6x can be directly related to the reduced matrix elements 
for electromagnetic transitions. Using well-known sum-rules for these matrix elements 
one finds a relation between the deformation length and the nuclear masses and sizes. 


For isoscalar excitations one gets [|I4 
52 = 27r 


5^2 = 


27r fr 


A (2A + 1) 


(42) 


niN < > AEx ’ 3 ttin '' ' ' AE^ 

where A is the atomic number, < > is the r.m.s. radius of the nucleus, and E^ is 

the excitation energy. 

The transition potentials for isoscalar exciations are 

dUopt{r) 


Uo{r) = 3Uopt{r) +r- 


dr 


(43) 


for monopole, and 

U.,{ r ) = , 

dr 

for quadrupole modes. 


(44) 
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For dipole isovector excitations 


5i = 


71 A 


(45) 


2 rriN NZ 

where Z (N) the charge (neutron) number. The transition potential in this case is 


14 


Ui{r) = X ( 


N 


A 


Z ^ ^ (IUqpi 


(Pa 


dr 3 dr^ 


opt 


(46) 


where the factor x depends on the difference between the proton and the neutron 
matter radii as 


X- 


2{N- Z) 


I {Rn + Rp) 


AR, 


np 


Rn 


(47) 


Thus, the strength of isovector excitations increases with the difference between the 
neutron and the proton matter radii. This difference is accentuated for neutron- 
rich nuclei and should be a good test for the quantity ARnp which enters the above 
equations. 

The time dependence of the matrix elements above can be obtained by mak¬ 
ing a Lorentz boost, assuming that U\ is the scalar part of a four-vector with zero 
vector-potential. Since the potentials Ux r(t) peak strongly at t = 0, we can safely 
approximate 6{t) ~ 6*(t = 0) = 7r/2 in the spherical harmonic of eq. (|4l|). One gets 

= 7 ^5= < > Kv(« = |) f/Alr(i)l . (48) 


where r{t) = 

Using the Wigner-Eckart theorem, the matrix element of the spherical harmonics 
becomes 

(2Jj -|- 1)(2A -|-1) 


< IfMf\Yxp\RMi >= 


47r(2Jj + 1) 



fIfXiA 

\—MfpMi) 

V 0 0 0/ 


(49) 


For high energy collisions, the optical potential U{r) can be constructed by using 


the t-pp approximation |^. One gets 


hv 


U{r) — —— axTN {o^NN + i) 


r') p2(r — rO “i 


(50) 
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where a^N is the nucleon-nucleon cross section, and a^N is the real-to-imaginary- 
ratio of the forward {6 = 0°) nucleon-nucleon scattering amplitude. A set of the 
experimental values of these quantities, useful for our purposes, is given in table 1. 

We are not interested here in diffraction and refraction effects in the scattering, but 
on the excitation probabilities for a given impact parameter. The strong absorption 
occuring in collisions with small impact parameters can be included. This can be 
done by using the eikonal approximation and the optical potential, given by eq. (^). 
The practical result is that the excitation probabilities for a given impact parameter 
b, including the sum of the nuclear and the Coulomb contributions to the excitation, 
are given by 


Pfi{b) = a^i{b) aj^(6) ^ exp { - aNN J J d^rpi{r') p 2 {r - r')} , (51) 


where r = ^/b‘^ + The corresponding excitation cross sections are obtained by an 
integration of the above equation over impact parameters. 


3 The effect of finite resonance widths 

Up to now we have assumed that the excited states are isolated states, with zero 
width. However, this assumption is not realistic and it is important to study the ef¬ 
fect of hnite resonance widths on the excitation amplitudes. This is specially relevant 
for the case of excitation of giant resonances, which have a broad structure. The sim¬ 
plest way to study this effect is by using the Coupled-Channels Born approximation. 
This approximation was used ref. to describe the excitation of the double giant 
resonance in relativistc heavy ion collisions. It is based on the idea that in such cases 
only the coupling between the ground state and the dominant giant dipole state has 
to be treated exactly. The reason is that the transitions to giant quadrupole and 
to the double-phonon states have low probability amplitudes, even for small impact 
parameters. However, an exact treatment of the back-and-forth transitions between 
the ground state and the giant dipole state is necessary. This leads to modihcations 
of the transitions amplitudes to the remaining resonances, which are populated by 
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the ground state and the GDR. In ref. 0 the application of the method was limited 
to the use of an schematic interaction, and the magnetic substates were neglected. 
These dehciencies are corrected here. The method allows the inclusion of the width 
of the giant resonances in a very simple and straightforward way. It will be useful 
for us to compare with the Coupled-Channels calculations with isolated states, as we 
described in the previous sections. Figure 1 represents our procedure. The CDR is 
coupled to the ground state while the remaining resonaces are fed by these two states 
according to hrst order perturbation theory. The coupling matrix elements involves 
the ground state and a set of doorway states >, where n specihes the kind 

of resonance and Ap are angular momentum quantum numbers. The amplitudes of 
these resonances in real continuum states are 


Q/0(g) =< 0(e) 


of > 


(52) 


where 0(e) denotes the wavefunction of one of the numerous states which are respon¬ 
sible for the broad structure of the resonance. In this equation e = — En, where 

E^ is the excitation energy and En is the centroid of the resonace considered. 

As we have stated above, in this approach we use the Coupled-Channels equations 
for the coupling between the ground state and the CDR. This results in the following 
Coupled-Channels equations: 


ih do{t) = > exp { - ^{Ei + e)t} a^^(t) 

= J de exp { - ^{Ei + e)t} , (53) 

and 

ih exp {i{Ei + e)t/h} ao(t) . (54) 

Above, (n = 1) stands for the CDR, oq denotes the occupation amplitude of the 
ground state and the occupation amplitude of a state located at an energy e 
away from the GDR centroid, and with magnetic quantum number p (p = —1, 0,1). 
We used the short hand notation V)(°^^(f) =< D^^\VEi,ij.{t)\0 >. 
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Integrating eq. and inserting the result in eq. (^), we get the integro- 
differential equation for the ground state occupation amplitude 

1 


ao{t) = -^2 / de |aW(e)p 

n 


X J dt' exp I - z(Ei + e)(t - «o(^0 5 


(55) 


where we used that = —oo) = 0. To carry out the integration over e, we should 

use an appropriate parametrization for the doorway amplitude Q!^^^(e). A convenient 
choice is the Breit-Wigner (BW) form 



Ti 

^ ^ 27r 

_e2 + rV4_ 


In 


this case, this integral will be the simple exponential 

I * |a"' (f) t exp { - } = exp { - 


. {El — iV/2)t 


h 


}■ 


( 66 ) 


(57) 


A better agreement with the experimental line shapes of the giant resonances is 
obtained by using a Lorentzian (L) parametrization for i.e., 


I««(e)t = - 

TT 


Ti E. 


(58) 


.(B|-Bf)2 + rlBh • 

where = Ei + e. The energy integral can still be performed exactly but now 
it leads to the more complicated result 

Jde |a*‘’(«)l^ exp{-i ^^‘^‘^* } = exp{-it^i—, 


(69) 

where AC{t) is a non-exponential correction to the decay. For the energies and widths 
involved in the excitation of giant resonances, this correction can be shown numer¬ 
ically to be negligible. It will therefore be ignored in our subsequent calculations. 
After integration over e, eq. (|55|) reduces to 


a„(() =-Si E y"'(() f_Ji [y<“>((')]’exp{ 


.(Bi-!ri/ 2 )(i-i') 
* s 


} ao(i') , 
(60) 
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where the factor iSi is iSi = 1 for BW-shape and iSi = 1 — iVi/2Ei for L-shape. 

We can take advantage of the exponential time-dependence in the integrand of the 
above eqnation, to rednce it to a set of second order differential eqnations. Introdncing 
the anxiliary amplitndes given by the relation 


ao{t) — 1 + X! ’ 


(61) 


with initial conditions A^(t = —cx)) = 0, and taking the derivative of eq. (|60|), we get 






Wt\t) 


,£1 

2 


A^{t) - 




l + Y.A,it)]. (62) 


Solving the above eqnation, we get ao(t). Using this amplitude and integrat¬ 
ing eq. (|54D , one can evalnate ae^i^(t). The probability density for the popnlation 
of a GDR continnnm state with energy in a collision with impact parameter b, 
Flip, Ex)I is obtained trongh the snmmation over the asymptotic {t —> cxd) contribn- 
tion from each magnetic snbstate. We get 


/ OO f p ^ 

dt' expjjB.i'} o„((') 


(63) 


where \a'^^\Ex — Ei)\^ is given by eq. (l5^) or by eq. (H), depending on the choice of 
the resonance shape. 

To hrst order, DGDR continnnm states can be popnlated throngh U2-conpling 
from the gronnd state or throngh Ul-conpling from GDR states. The probability 
density arising from the former is given by eq. (|63D, with the replacement of the line 
shape by its DGDR connterpart (defined in terms of parameters E 2 and 

r 2 ) and the nse of the appropriate conpling-matrix elements with the E2 

time dependence given by (29). On the other hand, the contribntion from the latter 
process is 


P2{h,Ex) = W^\Ex 

rt' 


/ OO 

jt' exp{iE,t'} {•£ (yy (('))' 


X 


df 


exp 


{ 


.iE^-iTi/2){t-t') 
* h 


} ao{t" 


, (64) 
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We should point out that eq. (^) is not equivalent to second-order perturbation 
theory. This would be true only in the limit ao(t) —> 1. In our approach, ao(t) ^ 1, 
since it is modihed by the time-dependent coupling to the GDR state. This coupling 
is treated exactly by means of the Coupled-Channels equations. We consider that 
this is the main effect on the calculation of the DCDR excitation probability. This 
approach is justihed due to the small excitation amplitude for the transition 1 —> 2, 
since ai{t) ao{t). 


Equations similar to can also be used to calculate the ISCQR and IVCQR 
excitation probabilities, with the proper choice of energies, widths, and transition 
potentials (e.g., VE 2 (t), or or both). 

In the next section we will apply the results of sections 2.1, 2.2 and 3, to analyse 
some examples of relativistic nuclear and Coulomb excitation. 


4 Applications 

We consider the excitation of giant resonances in ^osp]^ projectiles, incident on ^ospp 
targets at 640 A-MeV. This reaction has been recently studied at the CSI/SIS, Darm¬ 
stadt 1^. For this system the excitation probabilities of the isovector giant dipole 
(IVGD) at 13.5 MeV are large and, consequently, high order effects of channel cou¬ 
pling should be relevant. To assess the importance of these effects, we assume that the 
CDR state depletes 100% of the energy-weighted sum-rule and neglect the resonance 
width. The influence of resonance widths will be considered later, in section 4.2. 

4.1 Zero-width calculations 

As a hrst step, we study the time evolution of the excitation process, solving the 
Coupled-Channels equations for a reduced set of states. We consider only the ground 
state (g.s.) and the CDR. The excitation probability is then compared with that 
obtained with hrst order perturbation theory. This is done in hgure 2, where we 
plot the occupation probabilities of the g.s., |ao(t)P, and of the CDR, |ai(t)p, as 
functions of time, for a collision with impact parameter 6 = 15 fm. As discussed 
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earlier, the Coulomb interaction is strongly peaked around t = 0, with a width of the 
order At ~ b/jv. Accordingly, the amplitudes are rapidly varying in this time range. 
A comparison between the CC-calculation (solid line) and hrst order perturbation 
theory (dashed line) shows that the the high order processes contained in the former 
lead to an appreciable reduction of the GDR excitation probability. From this hgure 
we can also conclude that our numerical calculations can be restricted to the interval 
— 10 < r < 10, where r = {'jv/b) t is the time variable measured in natural units. 
Outside this range, the amplitudes reach assymptotic values. 

It is worthwhile to compare the predictions of hrst order perturbation theory 
with those of the harmonic oscillator model and the CC calculations. In addition 
to the GDR, we include the following multiphonon states: a double giant dipole 
state (2 (g) IVGD) at 27 MeV, a triple giant dipole state (3 g) IVGD) at 40.5 MeV, 
and a quadruple giant dipole state (4 g) IVGD) at 54 MeV. The coupling between 
the multiphonon states are determined by boson factors, as explained at the end 
of section 2.2. Direct excitations of the multiphonon states from the g.s. are not 
considered. The angular momentum addition rules for bosons yields the following 
angular momentum states: L = 0 and 2, for the 2 g) GDR state; L = 1, 2, and 3, for 
the 3 g) GDR state; and L = 0, 1, 2, 3, and 4, for the 4 g) GDR state. We assume 
that states with the same number of phonons are degenerate. In table 2, we show 
the resulting cross sections. The excitation probabilities and the cross section were 
calculated with the formalism of section 2. The integration over impact parameter 
was carried out in the interval bmin < b < oo. As we discuss below, the low-6 cut¬ 


off value bmin = 14.3 fm mocks up absorption effects. We have checked that 
the GG results are not signihcantly affected by the unknown phases of the transition 
matrix elements. Since the multiphonon spectrum is equally spaced, and the coupling 
matrix-elements are related through boson factors (see the end of section 2.2), the 
harmonic oscillator and the GG cross sections should be equal. In fact the numerical 
results of these calculations given in the table are very close. We also see that the 
excitation cross sections of triple- and quadruple-phonon states are much smaller than 
that for the 2 g GDR. Therefore, we shall concentrate our studies on the 2 g GDR, 
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neglecting other multiphonon states. 

Next, we include the remaining important giant resonances in Namely, the 

isoscalar giant quadrupole {ISGQ) at 10.9 MeV and the isovector giant quadrupole 
{IVGQ) at 22 MeV. Also in this case, we use 100 % of the energy-weighted sum 
rules to deduce the strength matrix elements. In table 3, we show the excitation 
probabilities in a grazing collision, with b = 14.3 fm. We see that hrst order pertur¬ 
bation theory yields a very large excitation probability for the IVGD state. This is 
strongly reduced in a c.c. calculation, as we have already disscussed in connection 
with the hgure 2. The excitations of the remaining states are also influenced. They 
are reduced due to the lowering of the occupation probabilities of the g.s. and of 
the IVGD state in the c.c. calculation. As expected, perturbation theory and c.c. 
calculations agree at large impact parameters, when the transition probabilities are 
small. For the excitation of the 2 ® IVGD state we used second-order perturbation 
theory to obtain the value in the second column. The presence ot the ISGQR and 
the IVGQR influence the c.c. probabilities for the excitation of the GDR and the 
2 ® IVGD, respectively. 

We should also consider the effects of strong absorption in grazing collisions, as 
discussed in section 2.2. In hgure 3 we plot the GDR excitation probability as a 
function of the impact parameter. In the solid line, we consider absorption according 


to eq. (|5l|) . In the construction of the optical potential we used the g.s. densities 
calculated from the droplet model of Myers and Swiatecki |^. As shown in ref. [|^, 
this parametrization yields the best agreement between experiment and theory. The 
dashed line does not include absorption. To simulate strong absorption at low impact 
parameters, we use = 15.1 fm as a lower limit in the impact parameter integration 
of eq. (^. This value was chosen such as to lead to the same cross section as that 
obtained from the solid line. 

In hgure 4, we plot the nuclear contributions to the excitation probability, and 
as a function of the impact parameter. We study the excitation of the isoscalar 
giant monopole resonance (ISGMR), the IVGDR, and the ISGQR. The ISGMR in 
^°®Pb is located at 13.8 MeV. As discussed previously, isovector excitations are sup- 
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pressed in nuclear excitation processes, due to the approximate charge independence 
of the nuclear interaction. We use the formalism of section 2.2, with the deformation 
parameters such that 100% of the sum rule is exhausted. This corresponds to the 
monopole amplitude ao = 0.054. The IVGDR and ISGQR deformation parameters 
are (5i = 0.31 fm and 62 = 0.625 fm, respectively. The IVGQR excitation probability 
is much smaller than the other excitation probabilities and is, therefore, not shown. 
The nuclear excitation is peaked at the grazing impact parameter and is only relevant 
within an impact parameter range of ~ 2 fm. Gomparing to hgure 3, we see that 
these excitation probabilities are orders of magnitude smaller than those for Goulomb 
excitation. Gonsequently, the corresponding cross sections are much smaller. We get 
14.8 mb for the isoscalar GDR, 2.3 mb for the ISGQR, and 2.3 mb for the IVGDR. 
The interference between the nuclear and the Goulomb excitation is also small and 
can be neglected. 

4.2 Effect of resonance widths 

We now turn to the influence of the giant resonance widths on the excitation dynamics. 
We use the GGBA formalism developed in section 3. Schematically, the GG problem 
is that represented figure 1. As we have seen above, the strongest coupling occurs 
between the g.s. and the GDR. 

In hgure 5, we show the excitation energy spectrum for the GDR, the DGDR (a 
short hand notation for the 2(DlVGD), ISGQR and IVGQR. The centroid energies 
and the widths of these resonances are listed in Table 4. The hgure show excitation 
spectra obtained with both Breit-Wigner (BW) and Lorentzian (L) line shapes. One 
observes that the BW and L spectra have similar strengths at the resonance maxima. 
However, the low energy parts (one or two widths below the centroid) of the spectra 
are more than one order of magnitude higher in the BW calculation. The reason for 
this behavior is that Goulomb excitation favors low energy transitions and the BW has 
a larger low energy tail as compared with the Lorentzian line shape. The contribution 
from the DGDR leads to a pronounced bump in the total energy spectrum. This bump 
depends on the relative strength of the DGDR with respect to the GDR. In hgure 6, 
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we show the ratio (TdgdrI(^gdr as a function of the bombarding energy. We observe 
that this ratio is roughly constant in the energy range Eiab/A = 200 — 1000 MeV 
and it falls beyond these limits. This range corresponds to the SIS-energies at the 
GSI-Darmstadt facility. 

We now study the influence of the resonance widths and shapes on the GDR and 
DGDR cross sections. This study is similar to that presented in ref. |P , except that we 
now have a realistic three dimensional treatment of the states and consider different 
line shapes. In the upper part hgure 7, denoted by (a), we show acDR as a function of 
T gdr^ treated as a free parameter. We note that the BW and L parametrizations lead 
to different trends. In the BW case the cross section grows with Tc^i/j while in the L 
case it decreases. The growing trend is also found in ref. @], which uses the BW line 
shape. The reason for this trend in the BW case is that an increase in the GDR width 
enhances the low energy tail of the line shape, picking up more contributions from 
the low energy transitions, favored in Goulomb excitation. On the other hand, an 
increase of the GDR width enhances the doorway amplitude to higher energies where 
Goulomb excitation is weaker. In hgure 7 (b) and (c), we study the dependence of 
(^GDR on Tgdr- In (b), the DGDR width is kept hxed at the value 5.7 MeV while 
in (c) it is kept proportional to <Jgdr-i fixing the ratio Togdr/^gdr = The 

first point to be noticed is that the BW results are sistematically higher than the L 
ones. This is a consequence of the different low energy tails of these functions, as 
discussed above. One notices also that <Jdgdr decreases with Tgdr both in the BW 
and L cases. This trend can be understood in terms of the uncertainty principle. If 
the GDR width is increased, its life-time is reduced. Since the DGDR is dominantly 
populated from the GDR, its short life-time leads to decay before the transition to 
the DGDR. 

To assess the sensitivity of the DGDR cross section on the strength of the matrix 
elements and on the energy position of the resonance, we present in table 5 the cross 
sections for the excitation of the GDR, DGDR, ISGQR and IVGQR, obtained with 
the GGBA approximation and 100% of the sum-rules for the respective modes. In this 
calculation we have included the strong absorption, as explained in section 2.2. For 
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comparison, the values inside parenthesis (and brackets) of the DGDR excitation cross 
section include a direct excitation of the L=2 DGDR state. We assumed that 20% of 
the E2 sum rule could be allocated for this excitation mode of the DGDR. The cross 
sections increase by less than 10% in this case. The value inside parenthesis (brackets) 
assume a positive (negative) sign of the matrix element for the direct excitation. 

Since the excitation of the DGDR is weak, it is very well described by eq. (|6^ and 
the DGDR population is approximately proportional to the squared strength of 1 /^ 2 ) 
Therefore, to increase the DGDR cross section by a factor of 2, it is necessary violate 
the relation of eq. (35) by the same factor. This would require a strongly anhar- 
monic Hamiltonian for the nuclear collective modes, which would not be supported 
by traditional nuclear models |^. Arguments supporting such anharmonicities have 
recently been presented in ref. . Another effect arising from anharmonicity would 
be the spin or isospin splitting of the DGDR. Since the Goulomb interaction favors 
lower energy excitations, it is clear that a decrease of the DGDR centroid would in¬ 
crease its cross section. A similar effect would occur if a strongly populated substate 
is splitted to lower energies. To study this point, we have varyied the energy of the 
DGDR centroid in the range 20 MeV < Edgdr < 27 MeV. The obtained DGDR 
cross sections (including direct excitations) are equal to 620 mb, 299 mb and 199 mb, 
for the centroid energies of 20 MeV, 24 MeV, and 27 MeV, respectively. Although 
systematics of the DGDR excitation do not show large deviation of the centroid 
energy, the data are not conclusive, and more experiments are clearly necessary. We 
conclude, that from the arguments analysed here, the magnitude of the DGDR cross 
section is more sensitive to the energy position of this state. The magnitude of the 
DGDR cross section would increase by a factor 2 if the energy position of the DGDR 
decreases by 20%, as found in ref. [|T^ . 


5 Discussion and conclusions 

In this paper we investigated at great length the excitation of giant resonances in 
heavy-ion reactions. Both the single- and double- giant dipole resonances were con- 
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sidered. The effect of the finite lifetimes of these resonances on their excitation 
probabilities was carefully assessed. The comparison with the available experimental 
data shows that some physics is still missing. Here we address this issue. 

In our discussion of the excitation of a damped giant resonance, the damping 
arises from the coupling to the large number of non-collective states that surround 
the GDR and shares with it its quantum numbers. One should keep in mind that 
our final result for the excitation probability involves an implicit average over the 
“chaotic” degrees of freedom whose quantum manifestition is just the fine structure 
states. At this point one is reminded of a well known fact in reaction theory, namely, 
ensemble or energy averaged cross sections contain two pieces: one obtained from 
an average amplitude, or “optical” piece, and a second piece which arises from the 
fluctuations. We expect similar contribution of the fluctuations to the excitation 
probability in the case of the GR. Here, however, the fluctuations are in the “host” 
nucleus and not in the compound nucleus. 

At this point, we recall similar type of fluctuations which constitute the dominant 


piece in the case of deep inelastic heavy ion reactions when it is assumed that 
only chaotic channels are involved in the inelastic transitions. The investigation of 
the effects of fluctuations on the excitation of giant resonances in heavy-ion reac¬ 
tions, following the procedure of |M, is underway and will be reported in a future 
publication. 


Acknowledgements 

We would like to express our gratitude to Dr. Hans Emling for useful comments 
and suggestions during the development of this work. One of us (GAB) acknowledges 
finantial support from the GSI-Darmstadt. This work was supported in part by the 
Brazilian agencies: GNPq and FINEP. This work was also partially supported by 
the National Science Foundation through a grant for the Institute for Theoretical 
Atomic and Molecular Physics at Harvard University and Smithsonian Astrophysical 
Observatory. 


24 











References 


[1] G. Baur and C.A. Bertulani, Phys. Lett. B174, 23 (1986); Phys. Reports 163, 
299 (1988); T. Aumann, C.A. Bertulani and K. Suemmerer, Phys. Rev. C51, 
416 (1995). 

[2] H. Emling, Prog. Part. Nucl. Phys. 3 (1994). 

[3] P. Chomaz and N. Frascaria, Phys. Reports 252, 275 (1995). 

[4] P.G. Hansen, Nucl. Phys. A553, 89c (1993); C.A. Bertulani, L.F. Canto and 
M.S. Hussein, Phys. Rep. 226, 282 (1993). 

[5] K. Alder and A. Winther, “Coulomb Excitation”, Academic Press, New York, 
1965. 

[6] L.F. Canto, A. Romanelli, M.S. Hussein and A.F.R. de Toledo Piza, Phys. Rev. 
Lett. 72, 2147 (1994). 

[7] C.A. Bertulani, L.F. Canto and M.S. Hussein, Phys. Letters B (in press). 

[8] A. Winther and K. Alder, Nucl. Phys. A319, 518 (1979). 

[9] J.D. Jackson, “Classical Electrodynamics”, (Wiley, NY), 1975. 

[10] A.R. Edmonds, “Angular Momentum in Quantum Mechanics”, Princeton Uni¬ 
versity Press, Princeton, New Jersey, 1960. 

[11] M.S. Hussein, M.P. Pato and A.F.R. de Toledo Piza, Phys. Rev. C51, 846 (1995). 

[12] J. Norbury and G. Baur, Phys. Rev. C48, 1915 (1993) 

[13] C.A. Bertulani and V. Zelevinsky, Phys. Rev. Lett. 71, 967 (1993); Nucl. Phys. 
A568, 931 (1994). 

[14] G.R. Satchler, Nuc. Phys. A472, 215 (1987). 

[15] M.S. Hussein, R. Rego and C.A. Bertulani, Phys. Reports 201, 279 (1991). 


25 



[16] L. Ray, Phys. Rev. C20, 1857 (1979). 

[17] W.D. Myers and W.J. Swiatecki, Ann. Phys. 55, 395 (1969); 84, 186 (1974). 

[18] C. Volpe, F. Catara, Ph. Chomaz, M.V. Andres and E.G. Lanza, GANIL preprint 
95/10, Nucl. Phys. A, in press; E.G. Lanza, private communication. 

[19] D.M. Brink, J. Neto and H.A. Weidenmuller, Phys. Lett. B80, 170 (1979). 


26 



Figure Captions 


Fig. 1 - Schematic representation of the excitation of Giant Resonaces, populated 
in heavy ion collisions. 

Fig. 2 - Time-dependence of the occupation probabilities |aoP and |aip, in a col¬ 
lision with impact parameter 6 = 15 fm. The time is measured in terms of the 
dimensionless variable r = {v'j/b) t. 

Fig. 3 - The GDR excitation probabilities as functions of the impact parameter, 
for sharp and smooth absorptions. 

Fig. 4 - Nuclear excitation probabilities as functions of the impact parameter. 

Fig. 5 - Excitation energy spectra of the main Giant Resonances for both Breit- 
Wigner and Lorentzian line shapes. 

Fig. 6 - Ratio between the DGDR and the GDR cross sections as a function of the 
bombarding energy. 

Fig. 7 - Dependence of (Jgdr and ancDR on the GDR width, treated as a free 
parameter. For details see the text. 
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Table 1. Parameters for the nucleon-nucleon amplitude, 
(fcATAr/dvr) aNN (* + «7VAr)- 


0 °) = 


E [MeV/nucl] 

a^N [fm^] 

OnN 

85 

6.1 

1 

94 

5.5 

1.07 

120 

4.5 

0.7 

200 

3.2 

0.6 

342.5 

2.84 

0.26 

425 

3.2 

0.36 

550 

3.62 

0.04 

650 

4.0 

-0.095 

800 

4.26 

-0.075 

1000 

4.32 

-0.275 

2200 

4.33 

-0.33 


Table 2 : Excitation cross sections (in milibarns) of the IVGDR, and of the n x 
GDR states in the reaction ^osp^ _|_208 g^g MgV.A. A comparison with Erst 

order perturbation theory and the harmonic oscillator is made. 


State 

1st pert. th. 

harm. osc. 

c.c. 

IVGD 

3891 

3235 

3210 

2 (g) IVGD 

388 

281 

280 

3 ® IVGD 

39.2 

27.3 

32.7 

4 ®IVGD 

4.2 

2.4 

3.2 


Table 3 : Transition probabilities at 6 = 14.3 fm, for the reaction ^osp^ p 208 p^ 
at 640 MeV.A. A comparison with Erst order perturbation theory is made. 


Trans. 

1st pert. th. 

c.c. 

g.s. —> g.s. 

— 

0.515 

g.s. —^ IVGD 

0.506 

0.279 

g.s. —ISGQ 

0.080 

0.064 

g.s. —^ IVGQ 

0.064 

0.049 

g.s. —^ 2 (g) IVGD 

0.128 

0.092 








































Table 4 : Centroid energies and widths of the main Giant Resonances in 



GDR 

DGDR 

ISGQR 

IVGQR 

Er (MeV) 

13.5 

27.0 

10.9 

20.2 

T (MeV) 

4.0 

5.7 

4.8 

5.5 


Table 5 : Cross sections in milibarns for the excitation of giant resonances in lead, 
for the reaction _l_208 p^ g^g ]\/[eY_A_ 


GDR 

DGDR ISGQR 

IVGQR 

2704 

184 (199) [198] 347 

186 
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